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We present new results for the renormalisation and subtraction constants for the four fermion AF — 2 
operators, computed non-perturbatively in the RI-MOM scheme (in the Landau gauge). From our preliminary 
analysis of the lattice data at f3 = 6.45, for the B° — B° mixing bag-parameter we obtain Bg GI — 1.46(7)(1). 
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Matrix elements of the four fermion (4f) op- 
erators play an important role in the studies 
of CP-violation in the Standard Model and be- 
yond. Of particular interest are the matrix el- 
ements describing the 5° - B a and K° - K° 
mixing amplitudes. Lattice QCD is at present 
the best suited method to compute such quan- 
tities non-perturbatively (NP). To keep the the- 
oretical uncertainties under control, it is of 
vital interest to renormalise the correspond- 
ing 4f-operators non-perturbatively. An addi- 
tional difficulty arises when working with Wil- 
son fermions on the lattice, namely all the par- 
ity even (AF = 2) operators mix among them- 
selves. Therefore, besides the overall renor- 
malisation, the spurious mixing should be sub- 
tracted. 

1. Non-perturbative renormalisation of 
four fermion operators 

To compute the values of the subtraction and 
renormalisation constants we use the method 
described in ref. |Q. Following those papers, 
we compute the Green functions (GF) of all 
AF = 2 operators, sandwiched by the quark 
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fields, in a specific gauge (in practice we opt for 
the Landau gauge). After projecting the am- 
putated GF's onto all independent Dirac struc- 
tures we get Tij (p 2 ), on which the following (RI- 
MOM) renormalisation condition is imposed 



rii(p 2 )U- 

,(0) 



= rgV) 



(i) 



where I^.' denote the tree level values of the 
GF's. We employed this prescription in ref. |Q 
to compute the renormalisation and subtraction 
constants for all four fermion operators. Here, 
we focus on the parity even ones, i.e., 



Oi 
5 



&7 M ? blul + ^7^759 ^759 
bj^qh/^q - bj^qbj^q 
bqbq - bj 5 qbj 5 q 
bq bq + bj 5 qbj 5 q 
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bo^qbonyq . 



(2) 



As mentioned above, the lattice regularisation 
with Wilson fermions breaks chirality, thus in- 
ducing the spurious mixing among all of the 
above operators. Therefore to renormalise the 
operator Oi, we should subtract the effects of 
this mixing, namely 

5 

d 1 ( t i) = Z( t ia) 1 (a) + ^A i (a)O i (a) . (3) 
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Figure 1. Z(p) as computed at four values of (3 and 
rescaled to ft = 6.45 as indicated in eq. (Q). The full 
line correspond to the perturbative running of the 
operator Oi (/Lt) at NLO in the RI-MOM scheme (see 
eq. (|) in which we use A^L = 250 MeV). 
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Figure 2. Subtraction constants A2-5 computed 
at {3 = 6.45 for several values of the renormalisa- 
tion scale fi. We observe that, to a good accuracy, 
9A2-5/cfyi ~ 0, as it should. 



The condition (Q) allows us to compute all sub- 
traction A2-5(a) and the renormalisation con- 
stant Z(fia). We perform such a calculation for 
4 different values of the lattice spacing, corre- 
sponding to /? = 6.0, 6.2, 6.4, 6.45. At each (3, 
we work with 4 values of the light quark mass, 
which allows us to extrapolate to the chiral limit 
(in which the RI-MOM scheme is defined). In 
this extrapolation a special attention is given 
to the subtractions of the pseudo-Goldstone bo- 
son (PGB) pole. These artefacts are the con- 
sequence of the fact that the operators in our 
Green functions are inserted at zero momentum. 
To subtract the effects of the PGB pole, we ex- 
tend the method of ref. Q , in a way described 
in ref. @. 

In fig. [l] we plot Z(fj,a), obtained for all of 
our values of [3 and then rescaled to (3 = 6.45 
according to, 



Z{p a') — R(a', a) Z(\i a) , 



(4) 



where R(a',a) is /i-independent if the discreti- 
sation effects of O(iJLa) are neglig ible. In fig. g 
we also show the renormalisation group running 
at NLO in perturbation theory, i.e., 



Z(n) = Z(hq) 
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q s (^o) - cts(fJ-) 
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where /3^ =0) = 11, 70 = 4 and J^/ =0) = 
8 log 2 - 1933/726 §. We have chosen ^ ~ 



3 GeV. We see that eq. (||) describes our data 
quite well. Moreover, the O(iia) effects ap- 
pear to be small, a feature already observed in 
the case of the renormalisation constants of the 
bilinear quark operators, computed with this 
same method ||. In fig. we show A2_5, com- 
puted at (3 = 6.45. We note that they are very 
weakly dependent on the renormalisation scale 
/x, as one expects in presence of small lattice 
artefacts. We also mention that their values de- 
crease with the lattice spacing a (i.e. with larger 
value of 0). After fitting each one of them to a 
constant, we obtain 

A 2 _ 5 = {-4.6(2), -1.2(1), 1.0(2), 0.4(l)}x 10- 2 . 

For a more detailed discussion, including com- 
plete list of numerical results, see ref. [||. 

2. B° — B° mixing parameter 

We now use the results of the previous section 
to compute B^ 1 (fi) which parametrises the ma- 
trix element of the operator 0\ (/z) as 



(B q \0 1 ( f i)\B q ) = -m 2 B J 2 Bq B x (/i) 



(6) 



Its value is extracted from the plateau of the 
ratio 
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Figure 3. Verification of the scale independence of 
the parameter B\ G1 defined in eq. (^). 
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Figure 4. Heavy quark (meson) extrapolation to 
1/mB,. The linear extrapolation is made with 4 
heavier mesons only, whereas for the quadratic one 
all 6 masses are used. 



where and P are the heavy-light axial cur- 
rent and the pseudoscalar density, respectively. 
Our sample contains 100 independent gauge 
field configuration generated on a 32 3 x 70 lattice 
in the quenched approximation. We work with 
4 light and 6 heavy quark masses. The directly 
accessed heavy-light pseudoscalar mesons (with 
light quarks around the strange quark mass) are 
in the range mp £ [1.7,3.6] GeV. In the above 
ratio we fix t x = 21 and find the plateau in 
the interval t y £ [45, 54], for all combinations of 
our light and heavy quark masses. We compute 
Bi{ijl) for 17 values of the renormalisation scale 
/i, which we then convert to the renormalisation 
group invariant form (RGI) at NLO accuracy, as 
To 



B* GI =a s (ri 2 Ml + 



47T 



J RI Bf J (M),(7) 



where /3q, 70 and Jpj are already defined after 
eq. <S). In fig. we show the quality of our 
B* G X 

In order to get the value of the B-parameter 
relevant for the B — B° mixing amplitude, we 
need to extrapolate to the physical point cor- 
responding to \/m,B d s ■ For that purpose we 
construct the quantity $(mp): 
7o 



$(mp)= 



a s (mp) \ 2/3 
a s (m B ) J 

mp 



a Q + 
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(8) 



where the argument in <fr(mp) indicates that 
the B-parameter is computed for the pseu- 



doscalar meson of mass nip. The role of 
the evolution factor in (^) is to resum the 
terms oc log (ms/mp), thus removing them 
from the expansion in l/mp. The extrapolated 
point, $(mp,), is simply B^ GI , the one that is 
needed for the B° — B° mixing. The 1/mp- 
extrapolation is illustrated in fig. |]. Our pre- 
liminary results read, 



Bif = 1.46(7)(1), 



B 



RGI 



= 1.38(3)(1), 



where the first error is statistical and the second 
is the difference between the results of the linear 
and quadratic extrapolations. 
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